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Abstract: The starting process with countable number of types fj,(t) generates a stopped 
branching process £(i). The starting process stops, by falling into the nonempty set S. It is 
assumed, that the starting process is subcritical, indecomposable and noncyclic. It is proved, 
that the extinction probability converges to the cyclic function with period 1. 
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1. Let us consider a measured state space (V, A), where A is the a-algebra on X. 
On this space we consider unbreakable homogenous Markov process with transition 
^ ■ probability P(t, x,A), where t denotes time, x £ X and A £ A. Considering every 

trajectory of the given process as an evolution of the motion of a particle, P(t, x, A) 
can be interpreted as a probability that a particle, which starts its motion from x £ X, 
falls into the set A £ A till the time t. It is assumed, that the time is discrete and the 
lifetime of a particle is equal to 1. At the end of its life the particle promptly gives 
rise to a number of offsprings, starting position of which are randomly distributed 
on the space X. The number and the position of these offsprings depends only on 



* Chair of Theoretical and Applied Statistics, Department of Mechanics and Mathematics, Ivan 
Franko National University of Lviv, Universytetska Str. 1, Lviv, 79000, Ukraine. 

tC.A.S.E. - Center for Applied Statistics and Economics, Ladislaus von Bortkiewicz Chair of 
Statistics of Humboldt- Universitat zu Berlin, Spandauer Strafie 1, D-10178 Berlin, Germany. Email: 
ostap.okhrin@wiwi.hu-berlin.de. 

* Chair of Theoretical and Applied Statistics, Department of Mechanics and Mathematics, Ivan 
Franko National University of Lviv, Universytetska Str. 1, Lviv, 79000, Ukraine. 



1 



the position of the particle- ancestor at the transformation time point. Further every 
offspring evolutes analogously and independently of other particles. 



Let Hxt(A) be a random measure, which for every A G A is equal to the number of the 
particles at time point t, types of which fall into set A, under condition that the process 
started with one particle x G X. fXt(A) is a random measure equal to the number of 
particles at the time t, which types are from the set A, but without any knowledge 
about starting group of particles. 

Further we assume, that the space X consists of a countable number of elements 
x 1: . . . , x n , This means that the set of types of particles {Ti, . . . , T n , . . .} is count- 
able. 

Based on the measure /i x t(A) we introduce the multivariate measure /-i xt {A) 



oo rii 

""k'ij 1 



. ,. , EE I^Xijt\ x m) j if %m G A 
»xt(A) = < i=oj=i 



0, else 



m=0 



where x = {x n , . . . , x lni ,x 2 i, ■ ■ ■ , x 2n2 , • • •}, G X is the j-th element i-th type. 

Let us denote jV = {0, 1, 2, . . .}, and respectively A/" °° is an infinite dimensional mea- 
sured space with elements x; t e A/o- 

Having P(t, x,A) let us introduce P(t,x, y), (x, y G A/" °°), where P is a probability 
that we obtain vector y till time t, assuming that we started from x. P could be 
rewritten in terms of /x xt as 

P(t,x,y) = P{ Mxt (X)=y}. 

Let £{i) = (5n, . . . , 5 in , . . .), where 5^ is the Kronecker symbol, <5 xy = f^i 4^,, ^(«) 
is the particle of the i-ih type. We also assume, that a b = a^a^ 2 • • - a^ n • • • , a! = 
ai!a 2 ! • • • a n ! • • • , a = a x H h a n H , af l1 = a^a* - 1) • • • (a, - 6, + 1). 

Definition 1. Functional 

W)) = *M = e«p{/m*),W 

zs called a generated functional of the random measure ji, where s(x) is a measured 
bounded function. 

Generated functional P(s) is always defined, when < \s(x)\ < 1 and integral f In s(x) /i(dx) 
exists. 



For our process the generated functional is given as 
h(t,s(-)) = Eexp j^lns(z)/z t (cfe) j, 



2 



where \x t is the same multivariate measure as fi xt but not taking into account any 
combination of the starting position of the process. Further we will consider the case 
s(-) = const = s = (si, S2, ■ ■ •)• K is easy to check whether the introduced generated 
functional is generated, as in the case of finite number of types (in that case it is not 
a functional but a function). 

Let us denote 



It is proved in [3], that the introduced generated function follows the main functional 
equation (V t, r = 0, 1,2,.. .) 

h(t + r, s) = h(t,h(r, s)). 

Let us fix the finite subset S C A/" °°, ^ S. Stopped or S -stopped multitype branching 
process is the process £ xi (X), defined for t — 1, 2, . . . and x e A/" °° by equations 



From this, for the S'-stopped process £ xt (X), points of the set S are additional states 
of absorption compared to the process /j, xt (X). The latter had only one point of 
absorption 0. In contrast to the process fj, xt (X), in the S'-stopped branching process 
£ xt (X) single particles in generation t multiplies independently following probability 
law defined by the generated functional h(-), only if £ xt (X) ^ S. If the random vector 
£ xt (X) falls into the set S, the evolution of the process stops. 

Since the process fx xt (X) is a Markov chain, then 



h\t,s) 
hP(t,s) 
h(t,s) 



((h^\t,s))^,(h^\t,s)f>,...), 
(h £ U(t,s),h £ W(t,s),...). 




AUX), ifW, 0<v<t, ^ V (X)<£S 

fi xu (X), if Vv, < v < u, H XV (X) i S, /x xu (X) G S, u < t. 



P(t 1 + t 2 ,a,(3)= P(ti,<x,7)P(t 2 ,7,P). 



For further needs, we also consider probabilities P(t, a, r), defined as 



P(t,a,r) 



P(l,a,r), t = l; 

^fc s P(l,a,(3)P(t-l,(3,r), t > 2. 



(1) 



It is easy to see, that P(l, a, r) is a conditional probability of the event 
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Let 

Qr(t) = P{Z nt (X) = r} 

be the probability of an extinction of the S-stopped branching process £ xt (X) into 
state r G S till time £, starting from state n G A/q 

2. Main facts. 

Theorem 1. For any n ^ S, n / 0, r 6 S, t > 1 ZioZds 
t 

?r n (*) = EE c «*W n ' a )' ( 2 ) 
cues i=i 

where coefficients cot r (t, Z) can be found from 

c ar (t+ 1,1 + 1) = car(M), (3) 

t-i 

c ar (t + l,l) = 5 ar -^P(i,a,r), (4) 

CQIr(l,l) = <W (5) 

Proof. Let 

r = min{t :/x nt (X) G S} 
be the moment of the first fall into S, then for t > Z 

P{^(X) = r, r = /} = P{^(X) = r} = P(Z, n, r). 
Applying ([T]) to P(Z, n, r), Z > 2, we get 

P(Z, n,r) = ^P(l,n,a)P(Z-l,«,r) 

Ot£S 

= P& n > a ^( l - 2, a, r) - ^ P(l, n, a)P(Z - 1, a, r). 

The first sum on the right hand side of this formula can be transformed similarly 

J2 P(2, n, a)P(l-2, a,r)=^ P(3, n, a)P(Z-3, a, r)-£ P(2, n, a)P(Z-2, a, r). 
a£s ol£s ct&s 

Making the same transformations in the sum Ylotgs P{h n ' Q: )P(Z — i, ct, r), we get 

i-i 

P(Z, n,r) = P(Z, n,r)- ^ ^ P(Z - z, n, a)P(i, a, r), (6) 

ctes i=i 

l = 2,...,t, P(Z,n,r)=P(l,n,r). (7) 
As g»(t) = Ya=i P{h n, r), from formulas ©,© we get □ 



Late on we will consider the process similarly to [1]. Let 

A 1 (x,D) = E{U(D)} 

be the first factorial moment, where £ x i(D) is such a random measure, which for each 
D G A is equal to the number of particles at time point 1, which types belong to set D, 
conditional on ^-stopped process. It also taken into account that at the beginning there 
was only one particle of the type x G X, what means € xi (D) = £xii(D). From 
the linearity of E we have Ai(x, D) = E{£ xl (D)} = A\(xi, D). It is important 
that D could be a vector or a set. 

Definition 2. Let Ai(x,D) = A(x,D) and 

A n+1 (x,D)= / A n (y,D)dA(x,y)= / A(y, D) dA n (x, y). 
Jx Jx 

It is assumed, that A (x,D) = 1, if x G D and A (x,D) = else. 

In [4] it is proved, that iterations of the operator A coincide with the first moments of 
£. This means, that for matrix of the linear operator A(t), with A,ij(i) = A t (xi,Xj), it 
holds that A(t) = A 1 will take place, where A = A(l). 

Let 

B t (x, D u D 2 ) = Eft^A) • U(D 2 ) - U(D! n D 2 )} 
be the second factorial moment. 

For further work we have to introduce some definitions, describing classes of branching 
processes (see [3]). 

Definition 3. Branching process in which all types form a single class of equivalent 
types is called indecomposable. All other processes are called decomposable. Branch- 
ing process is called fully indecomposable if the set of types could be split-up into two 
nonempty closed sets. 

Definition 4. An indecomposable discrete time branching process is called cyclic with 
period d, if the greatest divisor for all t, such that (A t (xi,Xi)) > 0, is equal to d. If 
d = 1 then the process is called noncyclic. 

Definition 5. An indecomposable discrete time branching process is called subcritical, if 
the largest eigenvalue (Perron's root) S of the matrix A is smaller than 1, supercritical, 
if 5 > 1 and critical if 5 = 1 and f(xi)B l j k i/(xj)v(xk) > 0, where B % - k is the matrix of 
the operator B, and f and v eigenfunction and invariant measure respectively which 
correspond do the Perron's root 5. 

Assumption 1. The kernel E£ xt (S) is assumed to be indecomposable, noncyclic and 
subcritical. 
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Correspondingly to the assumption [T] the operator A, which is defined by the kernel 
Fi{£ xt (D)} in the space of measurable functions and in the space of measures, has the 
eigenfunction /(•) and the invariant measure z/(-), such that 

„ oo 

/ f(y)A t (x,dy) = f(x) = y^f{y i )A t {x 1 y i ), 

J * *=1 

„ oo 

/ A t (x,Y)u(dx) = u(Y) = y2A t ( Xi ,Y)u(xi). 
Jx i=i 

Further we assume, that < X\ < f(x) < x% < oo, u(X) < oo and 

» oo 

/ f(y)u(dy) = l = y2f(y i )v(y i ). (8) 
Jx i=i 

The operator induced by the above defined kernel in the space of bounded functions 
has {1} as an isolated point of the spectrum. 

Assumption 2. We assume E{/j, £ ^ 1 (x i ) log ^s{i)ii x d} ^ s finite for Wi,j — 1,2, — 
Assumption 3. The expansion A t (x,y) = J2k f{ x k)^k v {yk) e2; 2sis. 

As in indecomposable, noncyclic, subcritical processes with discrete time all absolute 
values of eigenvalues are less than one, then based on the assumption [3] we can conclude, 
that when t — > oo 

A t ( Xi , yj ) = + o{Sl), 

where 5 is the largest eigenvalue. Thus 

lim A t (xi, y^S'* = f{x i )v(y j ). (9) 
t— yoo 

Let us denote 

R\t,s) = l-h*(t,8), 
R(t,s) = (R\t,s),...,R n (t,s),...), 
R(t,0) = Q(t) = (Q 1 (t),...,Q n (t),...) = \imR(t,s). 

s— >0 

As in the case with the finite number of types, the following inequalities could be easily 
proved (see [3]) 

< R\t, s) < Q\t) < \s\ < 1, (10) 
\R i {t,s)\<2Q i {t) 0<|s|<l. (11) 

(ITT1) implies that for the degenerating branching processes R l (t, s) converges uniformly 
to zero on < \s\ < 1. 

We need following technical assumption on the process 



6 



Assumption 4. Let A 1 > for some t > in the sense Vi,j a%j > and h l (t,s) ^ 

Aj(t). 

Hereafter the notation A = {aij} > 0, means that > Vz,j, and the notation 
A > B, where A = {ciij}, B = {%} are matrices, means that > Vi, j. 

Let h(s) = h(l,s). 

Assumption 5. Following the above defined assumptions for this process, it holds that 

l-h(s) = [A-E(s)](l-s), (12) 

where matrix E(s) with < s < s' < 1 satisfies conditions < E(s') < E(s) < A and 
\im s ^ 1 E(s) = 0. 

Theorem 2. With Assumptions EH2I 

r g!M ( \ 

urn — — — — = vvXi) 

t^oo f(x k )R k (t,s) V lJ 

uniformly on all s ^ 1, < s < 1. 

This theorem is proved analogically to theorem 1 on page 192 in [3], by replacing 
the right and left eigenvectors by eigenfunction and invariant measure respectively. 
Matrices are from the class of matrices of infinite measurable linear operator. 

Theorem 3. By assumptionsURMfor any i, j — 1, 2, . . . and for I — > oo probability that 
the process extinct to from one particle of type j over I is 

1- P(l,S(j),0) = K(S 3 )5 l (l + o(l)), K(S j )>0; (13) 

a) the limit of the conditional probabilities exists 

hmPUUX) = k|n^0}=p*, (14) 

t— >oo 

and the generating function h*(s) = ^ keA /-°° p^s^ is not depending on n and satisfies 
the relationships 

l-h*(h(-)) = 6(l-h*(s)), (15) 
/i*(0,...,0,...) = 0,/i*(l,...,l,...) = l; 

b) distribution p k has positive expectation 

h*(l) = Kmh*(s)= £ kjPL, 
keA/" °° 

where h*(s) = 

3 v ' as, 
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It is proved by mimicking the theorem 3 on page 198 from [3] with the use of theo- 
rem [2] for the representation of the limit of the generating function of the conditional 
distribution by getting result similar to one in [2]. 

Let us fix one more assumption 

Assumption 6. Let hij(s) = 9I £^ , then for all j, 1 < j < oo there exists such 

i, 1 < i < oo, that hij(0) are positive. 

From the equality 

MO) = P(0,£(i),£(j)) = P{^ (I)1 (I) = £(j)} 
this means, that the corresponding probabilities P(0,£(i),£(j)) are positive. 
To proceed further we need following lemma 

Lemma 1. Under the assumptions\T^ the limit of conditional probabilities is positive, 
for all i = 1, 2, . . . 

lim P{^ t {X) = £{i)\n ± 0} = p* £(i) > 0, 

Proof. The generating function h*(s) = J2kP*k sk * n Theorem [3] satisfies the equation 
(1151) . If we replace in this equation s by h(s), and repeat this replacement t times, we 
get the equality 

l-h*{h{t,s)) = 5 t {l-h*{s)), (16) 

where h(t, s) is t-th iteration of the function implied by the main differential equation. 
By differentiating ffT6|) with respect to Sj at s — 0, we obtain 

oo 

5>*(/^,0))M*,0) = ^(0) = SY^. (17) 
1=1 

As all coordinates of h(t,0) converge to 1, for t — > oo, then by the theorem [3] we can 
find such T and Ci, that h*(h(t, 0)) > C\ > for t > T. According to the assumption 
[6] this implies that for all 1 < j < oo we can found such i, that hij(t,0) > 0. For all 

ii, i 2 ,... ,it+i holds 

t 

(t,0)>l[h n h H+1 (0). 
1=1 

Thus (IT7|) implies 

t 

S %j) ^ C ^H h ^ °) > 0, VI < j < oo. 
i=i 

□ 
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Theorem 4. By the assumption[J\the limiting extinction probabilities g" = lim^oo q™(t), 
Vn ^ S, r G S, can be written in the series representation 

oo 

^ n = 5ZS c «r%^,a), (18) 
1=1 aes 

where c ar = lim^^ co,^, = $ctr ~ J27=i P( u ' ar )- 

Proof. Probabilities (t) increase with t and are bounded above by 1. Then the limit 
g„ = lirn^oo <2£(f) exists. 

We can pass to the limits on the left and on the right hand sides of the formula (J2J), 
when t — > oo, as for all a, r G £ holds that P(Z, a, r) < P(Z, a, r) and Chebyshev 
inequality and assumption |3] imply that 

{oo 
5>c«(£0'))>i 

oo 

< ^E{^(^(J))} 

oo oo 

= E«^^/(l + o(l)). 

i=l j=l 

This means that series J^P(i, ck, r) and ^P(/, o:,r) converge to each other. This 
implies (HE]). □ 



As in [2] let us consider the asymptotic behavior of for n — > oo. 

Theorem 5. Let assumptions [DEI are fulfilled and lim TT _ ! . 00 (r;,j/n) = Oj ; where a = 
(ai, a 2 , . . .). In this case for r G S and n — > oo 

g r n -Pf (log, n)^0, (19) 

where H(x) is a cyclic function with period 1, defined through the following equalities 

ro 

H{x) = $>P,(x), 

3=1 



L=—oo 



where (a, if) = YliLi a iKi, &i as ^ n (GIP; r o = max{r = r\ + r 2 + . . . : r G S}. 
Constants Cj = Cj(r, a, p*) depend on r, a and the limit distribution p* = {p* k } which is 
defined in lemmaUl 
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Proof. Let 9(1) = (9i(l), 92(1), ■ ■ ■) be a random vector, which components 9i(l) are 
equal to the number of particles of type i which give an offspring to the Z-th generation. 
Thus we can write, that for all a £ S, I > 1 n ^ S, we have 

P(l, n, a) = P i»nA X ) = «> d & = P) 

{f3:i<(3<et} 

P{9(0,l)=f3}P{ f i /3l (X) = <x\9(0,l) = (3}. (20) 

{f3-.i<(3<et} 
Under the assumptions of the theorem 

P{9(0,l) = f3} = ^^(P(l,£(t),0)) n ^(l-P(l,£( i ),0))^ 

= H ^^^ e -(a^(l+l+o(l)) (j + o(1) ) ) (21) 
/J. 

and the probability, not depending on n 
P{ti f3l (X) = a\9(0,l)=f3} 

oo 0k 

= E UU p {»% ) M) = <* m \m^o} 

{Oi.W}k=lj=l 

oo 0k 

— ► e nn^> ( 22 ) 

{Q;(jfc)} fe=l J=l 



where jPO are branching processes, whith the same distribution as fi £ ^ k ^ l (X). 

The summation in ^{cut^)} * s done over all such a^ k \ which YlkLi Sf=i — a - 
The statements in f[2"Uj) - ([2"2"j) imply, that the general component of the series (|18p for 
n — > oo, I — y oo can be written in the form 

ro — 

E ^K/3)E 5(WoS4lT) ^ 

Q:e5 {/3:i<^<a} 
xexp{ - (sl,K)8 Wo ^(1 + o(1))}, (23) 

where g(ot,f3) is an independent of n and I function. It is easy to see that in formula 
(fTSl) for n — > oo each component of series with any I > 1 converges to zero. 

Let us choose Li < L 2 in such way that sums 

oo L 2 

J2 S P L e-^ sL and £ S^e'^ 6 " (24) 

L=L2 L=— oo 

are small. We set Zj + log^ n = Lj + Xjn, for i = 1, 2, where < < 1. ( 12"3"|) and ( I2"lj) 
imply, that we can choose such L\, L 2 and no, that tails of the sum in f JT8|) . bounded 
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from 1 to l\ and from l 2 to infinity, are less then e/2, where e > is small. Elements of 
the series (j!8p with lx < I < l 2 can be replaced by a limited expressions (I2"5j) for n — > 00 
as well as for I — > 00. The number of summands in the sum X^ITh-i m expression f ll8p 
is finite l 2 — h — 1 = L 2 — L\ — 1. This means that no can be chosen in such a way, 
that for all n > no the approximation error will be also less than e/2. This implies the 
statement of the theorem, while e > is any real number. □ 



From the theorem it cannot be concluded directly, whether the coefficients Cj in the 
formula (1T9]) are such, that H(x) > 0. For this we introduce the next lemma. 

Lemma 2. Under assumptions\T^ there exists such a constant > ; that for some 
number uq 

g" > 6, for Vn with n > n and Vr E S. 



Proof. As for any t, g" = lim^oo g" (t) > g" (t), it is enough to prove, that the inequality 
Qr {t) > © > holds for any large enough t, for all r 6 S and n from n > n . Let 
us use the upper defined random vector 6(0, t) and introduce one more random vector 
6[{t — 1) = (9[(t — 1), 6' 2 (t — 1), . . .), where Q[{t — 1) is the number of starting particles 
of i-th type, with nonempty offspring set in the (t — l)-th generation, but empty in the 
t-th generation. For n = (ni,n 2 , . . .), ni > r + 1, where r = max r6 sr, we use the 
inequality 



Qr(t) = P{£ n4 (X) = r} 

> P{n at (X) = T,e'(t-l) = (r + l-T),e(0,t) = T£{l)}. (25) 

The right side of (|25|) we write product of Pi(n,t)P 2 (t), where Vi(n,t) = P{0'(t- 
1) = (r + 1 - r)£(l), 0(0, t) = rE(l)} and depends on n and t, but P 2 0O = 
P{n nt (X) = r\6'(t - 1) = (r + 1 - f)£(l), 0(0, t) = ff(l)} and depends only 
on t. From the definition of the random vectors 0(0, t) and 9'(t — 1) we have, that 



Vx(n,t) 



(ni - r - l)!(r + 



T _[p(t-l > £(l) J 0)] 



ni— ro-l 



(1 - P(t - 1, 5(1), 0)) ro+1 - ? (l - P(t - 1, 5(1), 0)) r 

00 

l[[P(l, £(1), 0r}P{fi' ro+1 _- rt (X) = I r + 1 - r ^ 0}; 



00 r fc 



k=i j=i 



(26) 
(27) 



Here /x, /x', /x( jfc ^ are branching processes, whose evolution is defined by a generating 
function h(s) = (^i(s), h 2 (s), . . .). Setting t — >■ 00, in such a way, that — > V > 
for n — >■ 00, we get in the right side of the equality (f2"6"|) a positive constant multiplied 
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by a conditional probability, which stays at the end of formula. Using the limiting 
relationship P{fi' t (X) = k | k ^ 0} — > p£, of the theorem [3] and the equality 

£ Pt(p ki a,£a),o)p k2 a,m,o)---) = h*(h(o)) 

keA/" °° 

we have that this conditional probability is equal in limit to h*(h(0)). Expression (1271) 
does not depend on n and is equal to the product n^itP^TS f° r t ~ * 00 • From 
lemma [1] this product is positive. That completes the proof. □ 
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